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' Abstract 

> : 

•i-H . 

^ ■ The results of an exploratory lattice study of heavy baryon spectroscopy are pre- 

■ sented. We have computed the full spectrum of the eight baryons containing a single 

heavy quark, on a 24^ X 48 lattice at (3 = 6.2, using an 0(a)-improved fermion action. 
We discuss the lattice baryon operators and give a method for isolating the contribu- 
tions of the spin doublets and to the correlation function of 
the relevant operator. We compare our results with the available experimental data 
and find good agreement in both the charm and the beauty sectors, despite the long 
extrapolation in the heavy quark mass needed in the latter case. We also predict the 
masses of several undiscovered baryons. We compute the A — pseudoscalar meson and 
S — A mass splittings. Our results, which have errors in the range 10 — 30%, are in 
good agreement with the experimental numbers. For the E* — E mass splitting, we 
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find results considerably smaller than the experimental values for both the charm and 
the beauty baryons, although in the latter case the experimental results are still pre- 
liminary. This is also the case for the lattice results for the hyperfine splitting for the 
heavy mesons. 



1 Introduction 



The discovery of the baryon at LEP |1[] and the claims of indirect evidence for A;, and 
Hfc semileptonic decays 0, have triggered an increased interest in the spectroscopy and 
weak decays of heavy baryons. The interest in the spectroscopy, in particular, has been 
considerably boosted after the announcement of the discovery of several spin-| charm and 
beauty baryons 0, ^ . 

The properties of hadrons containing a heavy quark can be studied using lattice QCD cal- 
culations, which provide non-perturbative, model-independent results. Experience gained 
through studies of heavy mesons has provided the framework for an investigation of the 
phenomenology of heavy baryons. Furthermore, the study of the spectrum of heavy baryons 
is a necessary precondition for the measurement of the weak matrix elements of semileptonic 
decays of beauty baryons. The computed masses and matrix elements can then be combined 
with an analysis carried out in Heavy Quark Effective Theory (HQET) |^ to extract an 
independent measurement of the CKM matrix elements Vcb and Vub- 

The subject of spectroscopy has been widely discussed in the literature, mainly using poten- 
tial models 0], HQET §], or a combination of the latter with chiral perturbation theory |p. 
Recently, there have been attempts to compute the mass of the Ah (one heavy quark and 
two light quarks) [^] and of the Ehh (two heavy quarks and one light quark) on the lattice 
TT| . In this paper, for the first time, the full spectrum of the lowest-lying baryons containing 



one heavy quark is computed. In particular, we define operators suitable for the simulation 
of baryon spin doublets with total spin ^ and |, like the (S,S*), (S',S*) and The 
quality of the signal we have observed and the agreement of our estimates with the available 
experimental data are good, thus giving us confidence in the reliability of our predictions. 

The quark content and quantum numbers of the baryons we have considered are summarised 
in Table |1|. On the lattices available at present, it is not possible to simulate directly the 
beauty quark, whose mass is larger than the cutoff. Therefore we have computed four heavy 
quark masses around that of the charm quark and interpolated (extrapolated) the results to 
the charm (beauty) quark, relying on the predictions of HQET. The masses of the charm and 
beauty quark were fixed from the masses of the D and B mesons, respectively. The results 
of the extrapolation in the heavy quark mass, could be checked carefully, given the relatively 
large sample of masses available. On the other hand, we have only used two different values 
of the light quark mass, thus limiting our ability to perform a detailed analysis of the chiral 
behaviour. Our results for the masses are given in Table |], where the first set of errors 
is purely statistical and the second set is an estimate of the systematic uncertainty in the 
calibration of the lattice spacing. Our results are in good agreement with the experimental 



determinations, where available^. We also present estimates for the A — pseudoscalar meson 
mass sphttings. Our resuhs are 

Ma, + Mb - 7 Ma, + Mb " ^ ^ ^ 

to be compared with the experimental values^ 

Ma -Md . , Ma, - Mb 

-i^ -^ = 0.100 3 -p -^= 0.033 5. 2 

Ma^ + Md ^ ' Ma, + Mb ^ ' ^ ' 

Similarly, for the S — A splitting, we find 

My: - Ma +9 Mg, - Ma, +5 

—±s ^ = 0.039 + I ^ = 0.017 + I 3 

Ms, + Ma, - 9 Ms, + Ma, " ^ ^ ^ 



which compares well with the experimental numbers 

My - Ma My - Ma 

The last number in Eq. (^), extracted from the data presented in ref. [|], is still preliminary. 

We also make a first attempt to estimate the spin splitting of the doublets (S*, S), (S*, S') 
and {Q*, Q) by isolating the contributions which the two particles give to the same correlation 
function. We find small, negative splittings, which, in most cases, become compatible with 
zero after the extrapolations because of the increased statistical errors. The simple quark 
model expectation is that the splittings are positive, although some of the experimental data 
are still inconclusive. If this expectation is confirmed by experiment, we could be facing a 
situation similar to that of the hyperfine splitting in heavy meson systems, where the splitting 
is underestimated using both the standard Wilson action and the Sheikholeslami-Wohlert 
(SW) action |]14|. The meson hyperfine splitting is sensitive to the chromomagnetic moment 
term which appears at 0(a) in improved fermion actions |]15[. We plan to investigate the 
sensitivity of the baryon hyperfine splitting by using the tadpole-improved Sheikholeslami- 
Wohlert action |jl6| . 



This paper is organised as follows: in Section 2 we discuss the baryonic operators which 
have been used in the present study and give details of the simulation. In Section 3 we 
explain our analysis procedures for the extraction of the masses. The measurement of the 
mass splittings is reported in Section 4. Our results and the comparison with the physical 
values are reported in Section 5. Finally, we present our conclusions. 



"'^The experimental evidence for the baryon is based on a collection of 11 events, and no estimate of 
the statistical error is given, see ref. [|l2j. We note that the physical S'^ and Eh states are mixtures of the 
states which we measure here, where the light-quark system has definite spin. It has been argued that such 
mixing (Tsj becomes negligible in the heavy quark limit |j. See the conclusions for further comments. 

^ Errors on the experimental data are added in quadrature. 
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Table 1: Summary of the quantum numbers of the eight haryons containing a single heavy 
quark. I, sj' are the isospin and the spin-parity of the light degrees of freedom and S, 
are the strangeness and the spin-parity of the baryon. 

2 Particles and Operators 

There are eight lowest-lying baryons containing one heavy and two light quarks (up, down or 
strange). The quantum numbers of the charm and beauty baryons are listed in Table |1], and 
their physical masses (see refs. p, ^, [1^, |18[), are given in Table 0. In the context of HQET 
at lowest order, it is possible to identify the spin-parity quantum numbers of the heavy quark 
and of the light system, within each baryon. Furthermore, heavy baryons with common light 
degrees of freedom exhibit common features; they are expected to be degenerate in mass, 
and to obey selection rules in weak decays. For example, the hyperfine mass splittings in the 
doublets (S*, S), (S*, S') and {Q*, Q) are expected to be 0(l/m/i), where nih is the mass of 
the heavy quark, and the weak semileptonic decay Sc is suppressed since it could only 

take place if the light quark system changed quantum numbers. 

2.1 Operators for Heavy Baryons 

The spectrum of the heavy baryons in Table [l| can be computed on the lattice by using the 
following three interpolating operators: 

a = eabc{r'^Cj,l")h''; O, = eabc{r^Cj,l'W; O'^ = eabc{r^Cj^l"')h'' (5) 

where a, b, c are colour indices, eabc is the totally antisymmetric Levi-Civita tensor, C is 
the charge conjugation matrix, /, /' are light quark fields, and h is the heavy-quark field. 
The (implicit) spinorial index of the three operators, is the (implicit) uncontracted Dirac 
index carried by the heavy quark field. By employing the operators in Eq. (^, one creates 
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Table 2: Heavy baryons considered in this project. Our results are quoted with a statistical 
error (first) and a systematic error (second) arising from the uncertainty in the calibration 
of the lattice spacing. Where available, we report the experimental data. 
("I" For the error on the mass, see footnote 1.) 

physical states whose heavy quark and light quark systems have definite quantum numbers, 
corresponding to the HQET description at lowest order. 

The operator corresponds to sf' = 0^ spin-parity for the light degrees of freedom and a 
total spin-parity for the baryon = The total isospin of the light degrees of freedom 
is / = if Z = M and I' = d, and / = | if one of the light quarks is the strange quark. 

Consider the two-point correlation function G^: 

G,ip, t) = Y: e-'^-%0,{x, t)O5(0, 0)). (6) 

X 

For large time separations, using antiperiodic boundary conditions in the time direction, this 
becomes in the continuum limit 

G,{p,t) g{e-^*(M + ^)-e-^(^-*)(M-^)} 

+ ^{e-''-^\Mpp-i)-e-^--^^-Hi + Mpp)} (7) 

where p^ = {E,p) is the 4-momentum of the baryon and p'^ = {E, —p). In Eq. (0), we have 
included the contribution of the parity partner (PP) baryon, with 4-momentum = {Epp, p) 
and q^ = {Epp, —p). This particle contributes to the correlation function because it has a 
non-zero overlap with the operator given in Eq. (^). At zero momentum, we choose 



an appropriate combination of spinorial indices such that the baryon, but not the parity 
partner, propagates forward in time. We find that the amphtude of the parity partner 
(Zpp) propagating backward in time is much smaller than that of the forward-propagating 
baryon {Z'^), and, in the following, we will neglect the contribution of the parity partners. 

The case of the operator (9^ is more involved than that of O5 since it transforms reducibly 
under the parity extended Lorentz group. It is the tensor product of a four vector and a Dirac 
spinor and thus transforms as (1, i) © (i, 1) © (i, 0) © (0, |) (in SU{2) © SU{2) notation). 
It can annihilate/create particles of spin-parity |^ and as well as these particles' parity 
partners. With the two-point function for the operator defined as 

Gf.u{p, t) = E e"*'"^ < 0^{x, t)a(0, 0) >, (8) 

X 

we find, in the continuum limit, for large values of t and using antiperiodic boundary condi- 
tions in time, 

-Z,Z,{i>,^^ + Mi/2)(Pif + ^2^1(1^1/2 - Mi/2)(P2f )m^(Pi/2)} 

+parity partners — antiparticles , (9) 

where pj = {Ej,p) and where parity partner contributions are obtained from the original 
particle contributions with the replacement Mj — > —Mpp^j while antiparticle contributions 
are obtained from the original particle and parity partner contributions with the replacement 
M — ^ — M, p —>■ —p and t ^ T — t in the exponent. For any given momentum p^, the 
quantities (P^/^)^!y(p) and {Pl-'^)^u{p)^ i,j = 1, 2, are the spin projection operators of ref. 
and are given by 

{P^^^)t,u{p) = g,,u - ^7m7z^ - ^ {"^IfiPu + Pt.lu'^) , 

[pU^),u{p) = ^P,P., (10) 
iPi2^)f.u{p) = -^^{p^.Pu-il^.Pv) 
{Pir)Ap) = -^^^{iPp^lv-PpPv) ■ 
They are orthonormal and idempotent 

(PlMPkiV = 5''5APi)l (11) 



where /, J take on values | or |. They satisfy 

7^(p3/2)/; = ; p^.iP^'^'' = = (P='/2)'^>. ; 

P.iPliy" = = {Pl/'rPu for z, J = 1, 2 (12) 
and have the following useful properties 

^{Pl/'r = MPl/'r^ +for^=j, -for^^j 

p((p3/2)M^^ _ (P3/2)M-J^. (13) 

The properties of Eqs. (|I2[) and (|I3|) guarantee that the spin-| contribution satisfies the 
appropriate Rarita-Schwinger equations and that the spin-| contributions satisfy the 
appropriate Dirac equations. The diagonal projectors are furthermore complete: 

9,. = (P'^XAp) + {PlhAp) + {PlhM ■ (14) 



To extract the masses of the spin-parity | and | particles, one needs only to compute the 
correlators (|^) at rest. In this case, the projectors P^/^ and Pi/^ simplify to 

(p3/2)., ^ ^^^-iyy; 2,J = 1,2,3, 

{Pli"r = ^7V, (15) 

and only act on the spatial components of G^jy(0, t), i.e. /i, v = 1, 2, 3. Since the components 
corresponding to the other projection operators do not contribute to the spatial components, 
Gij{6, t), it is clear from the properties of P^/^ and P^/^ given in Eqs. (|lT]) and (0) that the 
I"*" contribution can be isolated by considering {P^^^y^Gjk(0,t) and the 1^ contribution, by 

1 /2 ' ' ~* 

considering (P^^/ y^Gjk{0,t). The contributions of forward-propagating parity partners are 
suppressed by taking suitable combinations of spinorial indices as discussed after Eq. (^ for 
the case of 6*5(0, t). Those of the backward-propagating parity partners are naturally smaller 
because the time intervals over which they propagate are much longer for most values of t that 
we consider in analysing Gjk{0,t). Furthermore, both contributions are again empirically 
found to be suppressed by the fact that the overlaps of the parity partner states with the 
operator O'^ are orders of magnitude smaller and their masses slightly larger than those of 
the original particles. 

When space-time is approximated by a hypercubic lattice, full Euclidean 0(4) symmetry 
is reduced to symmetry under the hypercubic group. This reduction means that most ir- 
reducible representations of 0(4) and its covering group become reducible on the lattice. 
Fortunately, the representations which concern us here, (1, |), (|,1), (|)0) and (0,1), be- 
cause of their low dimensionality, remain irreducible on the lattice |3^. Furthermore, when 
restricted to the diagonal cubic subgroup (i.e. the lattice analogue of the rotation subgroup). 



(1, |) © (|, 1) decomposes into the reducible representation | © | while (|, 0) © (0, |) reduces 
to |, where | and | are themselves reductions to the cubic group of continuum spin-^ and 
spin- 1 representations. Thus, the space-time transformation properties of the operators O5 
and Ofj^ on the lattice are analogous to what they are in the continuum. Moreover, using the 
results of Ref. one can show that the cubic representations | and | mix only with the 
following continuum spins: 

^] .17 9 /3\ 3 5 7 ^^g^ 

2 /cubic 2 2 2 \2/ cubic 2 2 2 

Therefore, if one isolates correctly the cubic | and | contributions to G^^{p,t), one isolates 
unambiguously the contributions of the continuum spin-| and spin-| states in the large 
time limit (assuming, of course, that higher spin states are more massive). It should be 
emphasized that this isolation of the cubic representations must be done carefully because it 
is not known, a priori, which of the spin-| or the spin-| states is more massive. Fortunately, 
at zero momentum the continuum rest frame projectors given in Eq. (|T5|) are sufficient 
because they implement the Clebsch-Gordan decomposition of the product representation 
spin-l®spin-| into spin-| and spin-|, a decomposition which survives the reduction of SU{2) 
to the double valued cubic subgroup because irreducible representations of SU{2) with spin 
less than or equal to | are irreducible when reduced to that subgroup. It should further be 
noted that properties of operators and states under parity tranformations are unaffected by 
the discretization of space-time. 

A similar discussion applies to the operator O'^. In fact, the structure of the corresponding 
correlator in terms of quark propagators is the same as that of the only effect of the 
additional Wick contraction in the correlator is to change the overall normalization. 

Thus we have shown in the present section how to isolate the contributions of the different 
physical baryon states to the two-point functions of Eqs. (0) and In the following 

sections we shall use these procedures to determine the heavy-light baryon spectrum. To 
improve the overlap of the interpolating operators O^, and O'^ with the corresponding 
physical baryon states, we smear these operators as described in the Appendix. Though this 
smearing further complicates the isolation of the various physical baryon state contributions 
for non-zero two-point function momentum p (please refer to the Appendix for details), the 
baryon spectrum that we obtain here only requires an analysis of zero-momentum, smeared 
two-point functions to which the discussion of the present section applies unchanged. 



2.2 Details of the Simulation 



Our calculation is performed on 60 SU (3) gauge field configurations generated on a 24'^ x 48 
lattice sX j3 = 6.2, using the hybrid over-relaxed algorithm described in ref. |]22|. The quark 



propagators were computed using the 0(a)-improved Sheikholeslami-Wohlert action, which 
is related to the standard Wilson fermion action via 

= - E Qi^)F,Ax)a,M^), (17) 

where k is the hopping parameter. The use of the SW action reduces discretisation errors 
from 0{ma) to 0{asma) ^3, 24|, which is of particular importance in our study of heavy 



baryons, where the bare heavy quark masses are typically around one third to two thirds of 
the inverse lattice spacing. 

The gauge field configurations and light quark propagators were generated on the 64-node 
i860 Meiko Computing Surface at the University of Edinburgh. The heavy quark propagators 
were computed using the 256-node Cray T3D, also at Edinburgh. 



Statistical errors are calculated according to the bootstrap procedure described in |g2[, for 
which the quoted errors on all quantities correspond to 68% confidence limits of the distri- 
bution obtained from 1000 bootstrap samples. 

In order to convert our values for baryon masses and mass splittings into physical units we 
need an estimate of the inverse lattice spacing in GeV. In this study we take 

a-^ = 2.9 ±0.2 GeV, (18) 

thus deviating slightly from some of our earlier papers, where we quoted 2.7 GeV as the 
central value p2|, |2^, 0. The error in equation (|1^) is large enough to encompass all our 
estimates for a^^ from quantities such as mp, /^r, ttin, the string tension ^/K and the hadronic 
scale Rq discussed in [^. This change is partly motivated by a recent study using newly 



generated UKQCD data ||2^; using the quantity Rq we found a ^ = 2.95 GeV. Also, the 



non-perturbative determination of the renormalisation constant of the axial current yielded 
a value of Za = 1.05(1) which is larger by about 8% than the perturbative value which 
we used previously. Thus the scale estimated from f.^ decreases to around 3.1 GeV which 
enables us to reduce significantly the upper uncertainty on our final value of a'^ [GeV]. 

Light quark propagators were computed for quark masses around the strange quark mass, 
corresponding to hopping parameters k = 0.14144 and 0.14226. Because each heavy baryon 
contains two light quarks, we can form three baryon correlators for each heavy quark mass, 
of which two have degenerate light-quark masses and one has non-degenerate light-quark 
masses. The masses of the light pseudoscalar meson which are needed for this study, were 



obtained in ref. |2^. Results extrapolated to the chiral limit (corresponding to a hopping 
parameter Kent = 0.14315 ~^'^) and to the mass of the strange quark (ks = 0.1419 ^|) are 
also tabulated there. 



The heavy quark propagators have been computed for four values of the heavy quark mass 
around that of the charm quark, corresponding to Hh = 0.133, 0.129, 0.125 and 0.121. The 



masses of the heavy-light pseudoscalar mesons can be found in ref. . 

In order to enhance the signal for the baryon correlation functions, the light and heavy quark 



propagators have been computed using the Jacobi smearing method [§0] , either at the source 
only (SL) or at both the source and the sink (SS). Since smearing is not a Lorentz-invariant 
operation, it might alter some of the transformation properties of non-scalar observables. We 
have found that such an effect is evident in the baryonic correlators at non-zero momentum, 
and we present the results of a study of these effects in Appendix ^ This issue, which does 
not affect the spectrum, represents an important new effect which is crucial in the extraction 
of the amplitude Z, and therefore in the measurement of the weak matrix elements entering 
the semileptonic decays of the A^. 



3 Analysis Details 

It follows from Eq. (0), that, for t > 

[G'5(0,t)]n = [G5(0,t)]22 = -[G'5(0,T-t)]33 = -[G5(0,T-t)]44. (19) 

Therefore, we define the A correlation function as 

GaW = \ [[G'5(0,t)]n + G'5(0,t)]22 - [6-5(0, T-t)]33 - G'5(0,T-t)]44] - Zl e-^^' (20) 

Similarly, for the E and S*, we define the correlation functions by taking suitable combina- 
tions of the equivalent components, after projection with the operators given in Eq. ([151) . 

3.1 The Effective Masses 

In Fig. |I] we show effective mass plots of the A and S baryons, in both the SL and SS cases. 
We compute the effective masses assuming that the correlators' time evolution is given by 
an exponential: 

This is justified since we have checked that the contribution of the parity partners propa- 
gating backward in time is completely negligible. 

The effective mass is smoother for SL than for SS correlators, because the former are more 
correlated in time. To establish a fitting range, we fitted the correlators to a single expo- 
nential in the range [tmin, ^max], where tmax was fixed at 21 for SS and 23 for SL correlators, 

■^A is a conventional name, by which we mean the baryon whose hght quarks are in a = 0+ state. 
Depending on the flavour of the latter, this baryon is either the physical h.{ud) or the S(us) with spin for 
the light quarks. A similar convention is used for S and S*. 



and tmin was varied between 8 and 19. The fits at fixed values of the hght and heavy kappas 
were obtained by minimizing the computed using the full covariance matrix. 



Kh = 0.129 and kh = K12 = 0.14144, for both SS and SL correlators. The behaviour of the 
correlator for the S* baryon is similar, and the features described below are common to all 
the masses considered in this study. 

By considering the x^/dof of the fits, as well as the stability of the results under variation 
of tmin, we make the following observations: 

1. The masses obtained from the fits to the SS correlators are stable and the x^/dof are 
acceptable for tmin > H- For the SL correlators, the x^/dof are acceptable only for 
^min > 16. This behaviour supports the hypothesis that by smearing both the sink and 
the source one enhances the overlap with the ground state. 

2. The masses obtained from fits to the SL and SS correlators agree around tmin > 17 

3. As a general feature, we observe that the statistical errors increase with decreasing 
light quark mass. This effect is more pronounced for SL than for SS correlators. 

The conclusion is that there is good agreement between SS and SL data, even if, in the latter 
cases, the plateaux are shorter and the errors slightly larger. Thus we quote the results 
obtained by fitting SS correlators in t G [12, 21], using those obtained with SL correlators as 
a consistency check. 

3.2 Mass Extrapolations 

We obtain the masses of the eight charm and beauty baryons by extrapolating first in the 
light quark masses and then in the heavy quark mass. 

3.2.1 Extrapolation in the Light Quark Mass 

In order to perform the extrapolation to the chiral limit, we use the three baryon masses 
obtained from both degenerate (i.e. kh = k,i2 = 0.14144, or 0.14226) and non-degenerate 
(i.e. Kii = 0.14144 and K12 = 0.14226) light quark correlation functions. We assume that, in 
the chiral regime, Mbaryon depends linearly on the sum of the two light quark masses, that 
is, 



As an exampli 



.e, we show in Fig. |^ the results of this analysis for both the A and the S, with 




(22) 



with = (k;^^ + 1^12^) /2. This is supported by our data for the masses in the A and S 
channels, as shown in Fig. ^, and by previous studies of hght meson masses ^ . In the A 
channel, extrapolating to nn = K12 = Kent gives the mass of the Ah, while extrapolating nn 
to Kcrit and interpolating k,i2 to Ks gives the mass of the H/^. Similarly, from the S channel 
we extract the masses of the S^, the and the flh- 

The results of this analysis, obtained from SS correlation functions, are summarised in Ta- 
ble ^. By performing the same type of analysis on the SL correlators, we obtained essentially 
indistinguishable results. We note that the difference in the statistical errors of the SS and 
SL masses, already present in the fits at fixed k, is amplified by the extrapolation to the 
chiral limit. This confirms our earlier conclusion that by using SS data one obtaines more 
precise results. 



3.3 Heavy Quark Extrapolation 

The physical masses of the charmed and beauty baryons are obtained by extrapolating the 
four sets of data, computed at = 0.133,0.129,0.125 and 0.121. In performing these 
extrapolations, we have been guided by the HQET and have expressed the dependence of 
the baryon mass Mbaryon(^p) on the heavy-light pseudoscalar meson mass, Mp, through the 
following function 

Mbaryon(Mp) = Mp + C + ^ (23) 

IVl p 

where the two constants C and A are the parameters of the fit. The masses of the charm and 
beauty baryons are obtained for Mp = Md or Mp = Mb respectively. In Table ^ we report 
the results, corresponding to the SS case, in physical units. The numbers corresponding to 
the charm and beauty masses have been obtained assuming = 2.9 GeV. The quoted 
systematic error arises solely from the uncertainty in the scale and has been estimated by 
letting vary by one standard deviation about its central value. 

In these fits, the coefficient A, which quantifies the size of the l/ruh corrections, is of the 
expected size, i.e. A = 0(Aqqj-,), ranging from about (350MeV)^ to (540MeV)^, depending 
on the particular baryon and on the fiavour of the light degrees of freedom. Of course, the 
0{l/mh) corrections play an important role in the case of the mass splittings, (see section 4), 
while they contribute much less to the value of each mass. As a further confirmation of this, 
we have also set A to zero and verified that the results of the extrapolation are essentially 
indistinguishable from those presented in Table ^ although the are significantly higher. 
Finally we have used a function of the kind 

Mbaryon(Mp) = C + A' {Mp - Mq) (24) 

where a different slope in Mp is allowed, and we have obtained A' compatible with 1 in most 
cases. All of this confirms that heavy quark symmetry is very well satisfied here. Moreover 
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Table 3: Masses of the A, S and S* in lattice units obtained by fitting the SS correlators in 
t G [12,21]. Also shown are the corresponding masses after extrapolation in the light-quark 
masses, using Eq. 



the insensitivity of the results to the different modeHmg functions gives us confidence, not 
only in the interpolation to the charm mass, but also in the long extrapolation to the beauty 
mass. We stress, once more, the total agreement between SS and SL results, both for the 
final numbers and the features of the extrapolations. Two examples of the extrapolation to 
the heavy scale, corresponding to the A and the S are shown in Figs. ^ 



4 Calculation of mass splittings 

Once the value of Hh is fixed to correspond to the physical quark mass by matching Mp 
to either the D- or the S-meson mass, no large uncertainties are expected to occur in the 
measurement of other charm or beauty hadron masses, which are largely determined by that 
of the heavy quark. On the other hand, splittings in the masses arise from the dynamics of 
the light quarks and their interactions with the heavy quark. Their study provides a test of 
HQET as well as important information on the size of various systematic effects. 

The mass splittings are small quantities in comparison to the baryon masses themselves. 
Thus, they are affected by relatively larger statistical errors, as well as being more sensitive to 
the fitting and extrapolation procedures adopted. Therefore, a particularly careful analysis 
is required. Once more, we will quote results obtained by fitting SS correlators; the SL 
correlators give consistent results, although the statistical errors are appreciably larger. 



4.1 A— Pseudoscalar Meson Mass Splitting 

The A — pseudoscalar meson mass splitting is very precisely measured experimentally, ex- 
pecially in the charm sector. Therefore, as we use the pseudoscalar mass to fix the value 
of the heavy quark k, the agreement of the lattice value of M\ — Mp with experiment re- 
flects the extent to which our calculation properly incorporates the dynamics of the light 
degrees of freedom. The amount of computational effort devoted to this calculation, both 



with static [^-[Q and propagating Wilson fermions is testimony to its importance. 



We summarize the results obtained so far |32|, Q and compare them with the numbers from 
this study in Fig. ||. 

For the analysis of this splitting we need the correlation function of the pseudoscalar meson, 
which was determined in an earlier simulation using the same heavy quarks as this study. 



but with one additional light quark, corresponding to ni = 0.14262. We find that there is 

very little to be gained by determining the difference of the masses from the time evolution 

of the ratio of the A and the pseudoscalar meson correlator^, as opposed to obtaining it 

^Since the behaviour of the baryon and meson time shoe correlators is different close to the centre of the 
lattice, the ratio method is only safe if one excludes the last few timeslices from the fitting range. 
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Table 4: A-pseudoscalar meson mass splitting, in lattice units, obtained from the difference 
of the fitted masses. 



by subtracting the two masses determined separately. One reason is that the statistical 
errors associated with our measurement of the A mass are much larger than those of the 
pseudoscalar, and would dominate the final uncertainty on the splitting in any procedure. 
Moreover, when the pseudoscalar correlators were computed, only the heavy quark propaga- 
tors were smeared. Here the light quark propagators are also smeared, and consequently the 
pseudoscalar and baryon correlation functions suffer from different statistical fluctuations. 
Therefore, we measure the splitting by taking the difference of the baryon mass, obtained as 
described in the previous section, and the meson mass, fitted in the range t G [12, 21], as in 
ref . [0] . The results, for all the kappa values, are reported in Table |. 

We perform a chiral extrapolation of the mass differences at each heavy k value. Although we 
have simulated only two values of the light quark mass we have computed baryon correlators 
corresponding to two degenerate and one non-degenerate case; this last set of data, however, 
cannot be matched with the mesonic data and cannot be used in the chiral extrapolation. 
Hence, the chiral extrapolation is modelled by a linear function of the two degenerate light- 
quark points. Both our results for the A extrapolation (see section 3) and the evidence 
reported in ref. ||25| for the meson, justify this procedure. 

We performed the extrapolation to the physical pseudoscalar meson masses following two 
different procedures, in order to have a consistency check ( see figure |^). 



A The splittings are first extrapolated in the inverse heavy quark mass, according to the 
formula 

B C 
Ma-Mp]M = A+— — - + — 2--, 

J Mp{Kh) Mf,{Kh) 

keeping the light quark mass fixed. The linear and quadratic extrapolations produce 
indistinguishable numbers. Then the two values of this splitting corresponding to the 
two degenerate light-quark configurations are extrapolated to the chiral limit for each 
heavy- quark mass. 



B We employ the reverse procedure in which the mass splittings are first extrapolated to 
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experiment [MeV] 
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Table 5: Results for the A — pseudoscalar splitting, at the physical masses, corresponding to 
= 2.9 GeV. The two methods illustrated in the text produce essentially identical results, 
in excellent agreement with the experimental number, given in the last column. 

the chiral hmit, keeping the heavy quark mass fixed. The results from the subsequent 
linear and quadratic extrapolations in the heavy quark mass are once again compatible. 

We conclude that the behaviour of the A-pseudoscalar splitting is well represented by a 
linear function of both the light quark mass and the inverse heavy quark mass. The results 
in physical units that we quote in Table || are obtained under this assumption. 

4.2 E - A Mass Splitting 

We have obtained the S — A mass splitting for the various k combinations, both by taking 
the difference of the fitted masses and by fitting the time evolution of the ratio of S and 
A correlators. The numbers obtained with the two methods are in good agreement, but 
the second procedure yields appreciably smaller errors and is smoother in the chiral limit, 
proving that it is particularly appropriate when one compares two correlators with a similar 
structure. The results at each value of the computed masses and after the chiral extrapolation 
are given in Table |[ 

The dependence of the splitting on the heavy quark mass is extremely weak, suggesting that 
the l/rrih corrections to the masses of the two baryons must be very similar and nearly cancel 
in the difference. This feature makes it particularly simple to perform the extrapolation to 
the physical masses, as the fits to linear and quadratic functions of the inverse pseudoscalar 
meson mass are essentially indistinguishable. This is clearly visible in Fig. |^. We note that 
our results, presented in Table |^, compare very well with experiment. 

4.3 Spin Splitting 

In the HQET, the mass difference within the spin doublets (S,S*), (S',S*) and (i7,fi*) is 
due to the coupling of the chromomagnetic moment of the heavy quark to the light degrees of 
freedom. It is therefore suppressed by inverse powers of the heavy-quark mass and vanishes 
in the infinite mass limit. 
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Table 6: Estimates of M^, — Ma in lattice units, for the various kappa combinations, obtained 
with the ratio method. The extrapolations are linear, as they were for the masses themselves. 



Because the splitting is such a small number, it is difficult to measure from our data using 
either the ratio of the correlators or the difference of the ffited masses. We find that the 
difference of the masses gives smoother chiral behaviour, and so we obtain our estimates 
using this method. 

We present our measurements of the splittings for each k combination in Table || and the 
extrapolated values at the physical masses in Table |^. They are negative within two standard 
deviations at fixed light-quark mass, but being affected by large statistical errors, become 
compatible with zero in the chiral limit. We obtain the splittings at the physical masses, 
extrapolating in the inverse heavy-quark mass, according to the following function 

Ms* - Ms] (Kh) = — ^ + jj^. (25) 

J Mp[Kh) Mf,{Kh) 

The splitting has been constrained to vanish in the infinite heavy-quark mass limit, as 
predicted by the HQET. The two extrapolations obtained either by including the quadratic 
term or by setting B = gave essentially indistinguishable results, as can be seen from Fig. |^, 
and good x^- We also checked the consistency of our data with the predicted behaviour, by 
adding a constant to the function ([25|), i.e. by allowing the spin-splitting to have a non-zero 
intercept, at 1/Mp = 0. We find that the values of the intercept are always compatible with 
zero, being of the order -30 MeV, with errors of about one hundred. 

Once more, the results were perfectly consistent. The values presented in Table |^ correspond 
to Eq. (p5D, with B = since this parametrization fits the lattice data very well and makes 
full use of heavy-quark scaling relations. 

We note that our value for the S* — Sc splitting is inconsistent with experiment, and in most 
cases our measured values are consistent with zero. This feature resembles the well-known 



puzzle of the spin splitting in the mesonic sector , whose resolution is believed to lie in a 



combination of discretisation and quenching effects, as argued in ref. |15]. In the light of this 
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Table 7: Baryon-pseudoscalar meson and baryon-haryon mass splittings in MeV. The avail- 
able experimental data are also shown, together with the corresponding references. The ex- 
perimental errors on the — Sc and S* — splittings are not published. 

and also considering the very large statistical errors, firm conclusion about the consistency 
of our results with heavy-quark scaling laws cannot be drawn. 

5 Physical Results 

In this section, we present a summary of the results obtained in this study, in a form which 
is easily comparable with the experimental data. All masses are given with an asymmetric 
statistical error arising from the bootstrap analysis, and a systematic error due solely to the 
uncertainty in the scale (see Eq. ([TsD). 

In Table ^ we quote the charm and beauty baryon masses, together with the experimental 
values, where available. Our results agree well with the experimental data for the charm 
sector, and also for A;, and S;,, despite the long extrapolation in the heavy mass scale needed 
in these cases. This gives us confidence in the reliability of our predictions for the masses of 
the undiscovered charm and beauty baryons. The quality of the results at the beauty mass 
was certainly enhanced by the number of heavy quark masses available for this investigation, 
which allowed us to try different extrapolation procedures and to perform consistency checks. 
On the other hand, we only have a limited sample of light quark masses. Although the light 
extrapolations were always smooth and reasonable, the chiral behaviour should be confirmed 
by using a larger number of light quark masses. 

We present the mass splittings of Section 4 in Table |^. In those cases where comparison 
with experiment is possible, we also compute the ratios of the splitting to the sum of the 
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Table 8: Estimates of Ms* — in lattice units at the various masses, obtained by taking 
the difference of the fitted masses, in lattice units. The extrapolations are linear, following 
the same procedure adopted for the other splittings. 

masses, to eliminate most of the uncertainty in the scale. These results were presented in the 
Introduction, see Eqns (|I|-^). The residual systematic uncertainty, which is always smaller 
than the statistical error, was not quoted. 

In those cases where a meaningful comparison with experiment is possible the agreement 
is very encouraging. Unfortunately, the mass differences, being small, are affected by large 
relative errors varying between 10 and 30%. Nevertheless, we stress the beautiful agreement 
with the experimental data, both at the charm and at the beauty mass. In particular, in 
our calculation of the Ah — pseudoscalar meson splitting, the agreement with experiment 
has significantly improved on previous calculations, performed with the standard Wilson 
action. We believe that this success is further evidence of the advantages of using the 
Sheikholeslami-Wohlert fermion action. 

6 Conclusions 

In this paper we have presented the result of a lattice study of heavy baryon spectroscopy. 
The spectrum of the eight lowest-lying heavy baryons, containing a single heavy quark, can 
be computed using the three baryonic operators in Eq. (|^). In addition to the calculation of 
the A and S masses we have discussed how to compute the spectrum of the spin doublets, 
(ZI,S*), (S',S*) and (f2,fi*), by isolating their contribution to the correlation functions of 
the operators and O'^. 

The computation of the mass spectrum proved feasible; the operators we have used have a 
good overlap with the various baryon ground states, in part thanks to the smearing both at 



the source and at the sink. Moreover, the extrapolations in both the heavy and hght quark 
masses are always smooth. The agreement between our estimates of the baryon masses and 
the experimental values is good, in both the charm and the beauty sectors. 

The computed A— pseudoscalar meson mass splitting is in good agreement with experiment, 
in contrast with the results of previous calculations performed with the Wilson fermion 
action. We believe that this is largely due to the use of the 0(a)-improved action to remove 
systematic effects. A similar positive conclusion can be drawn for the E — A splitting, 
although the statistical errors are still of the order of 25 — 30%. 

Our results are also in agreement with the predictions obtained with other non-perturbative 
methods both for the masses themselves and for the S — A splitting. In the case of 

the S' — S splitting, which is of the same nature as the E — A splitting, the calculation is 
complicated by the mixing arising between the two particles, whose total quantum numbers 
are the same. It has been noted that such a mixing, negligible in the heavy quark limit, 
would have the effect of increasing the splitting. Both our prediction and that of Savage 
are higher than experiment |T^. The disagreement would hence get even worse if we were 
to take the mixing into account. We stress, anyway, that this experimental result is still to 
be confirmed. 

In this exploratory study the masses have been determined with reasonable precison, but 
further studies are required to reduce both the statistical and systematic errors. The results 
presented in this paper are very encouraging and it looks likely that it will also be possible 
to measure the baryonic matrix elements. 
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A Effect of smearing on baryon correlation functions 



In this appendix we propose a description of the effect of a non-Lorentz invariant smearing on 
a non-scalar operator. Interpolating operators are smeared to improve their overlap with the 
physical states one wishes to create or annihilate. We discuss the case of 2-point correlation 
functions, restricting ourselves to the case of spinorial operators which have overlap with 
a single type of spin 1/2 particle, like O^lx) defined in Eqn. (^. Numerical evidence for 
this effect is also presented. It should be noted that the breaking of the Lorentz symmetry 
manifests itself only when one considers correlators at finite momentum, and has therefore 
no relevance to the determination of the spectrum. 

Let us write the general expression for a local baryonic operator, 

Jp{x) = (^(x)rV'(x))^p(x), (26) 

and consider the case where it has overlap with spin 1/2 states, like O^. It can destroy one 
such state, according to the relation 

{Q\J,{OMr) = Zu^;\p). (27) 

In Eqs. (^,^), p and r are the spinorial and polarization indices, respectively, the antisym- 
metric sum over colour is understood and F is a suitable combination of gamma and charge 
conjugation matrices. Finally, the amplitude Z is a Lorentz scalar. 

In general, a smeared baryonic operator can be written as 

= E /(|y-x|)/(|^-x|)/(|^-x|)(^(y,t)rv^(^,t))v^p(^,t). (28) 

y,z,w 

Because the smearing is performed only in the spatial directions, Lorentz symmetry is lost 
and only spatial translations, rotations, parity and time reversal survive. Therefore, the 
overlap of the operator Jp{x,t) with the state \p,r) is given by the more general expression 

mj;mp,r) = [z,mu^^\p} + z,{\p\hou'^^\p}]^ (29) 

where the amplitudes Zi and Z2 may depend on the magnitude of the three-momentum of 
the state \p, r) , in accord with the restricted symmetries of the system. 

Let us consider the case of a SS two-point correlator, for large t 
G%{t,p) = Y.mj;{^,t)J:{6,0)]\0)e-'f^-' 

X 

m 



EE;^(0|^;(0,0)|g,r)(g,r|J:(0,0)|0)e-^-V^^^^^ (30) 

\g,r) X ^" ' 

E (zii\p\) + Z2{\p\hoy'\p)^^''\p)(Zim) + Z2{\p\ho] 



E{p) 



pa 



This expression can be conveniently rewritten as 



E + m — a'^{E — m) E + m + a'^{E — m) 



AE 



-I 



AE 



2a ^ ^ 



J pa 

(31) 

where Zs = Zi + Z2, a = [Zi — Z2)j{Z\ + Z2) and full Lorentz symmety is recovered when 
a = 1. Eqn. (|3TD exhibits the following features: 

• the exponential fall-off is not altered by the smearing. This was expected since a 
smearing function which only extends in the spatial directions preserves the form of 
the transfer matrix; 

• it has the correct limit for a 1 : (p, t) oc G^"'^ (p, t) ; 

• the degeneracy among the amplitudes of different spinorial components of the correlator 
is lifted; 

• all the terms proportional to a vanish for p — > 0, so that the effect disappears at zero 
moment umQ. 

In the following, we will check this effect against the numerical data. It is convenient to 
rewrite Eq. ( pl| ) in the form 

Zl 



~E{p)t 

~2E 



{mZjn + IqEZe - p- iZp] , with 



E + m + a^i^E — m) 



E + m — a'^{E — m) 



aZl 



2E 



s 

2m 



(32) 

We have found that Ze and Zm are compatible, as they should be, considering that they 
differ by terms proportional to / Em, which are very small for the values of momenta and 
masses in our simulation. Furthermore Zp is significantly different from Ze and Zm, which 
shows that a is different from one. Zg and a are given by 



EZe + mZf^ 



and 



a 



zr 



(33) 



^ E + m 

The results of this exercise are presented in Table ^ for the A baryon with momentum 
p = (^,0,0) and (^, ^,0), for masses corresponding to the four values of and kh = 
K12 = 0.14144. Using the estimated values of Z'^ and a we have also recomputed Ze and 
Zm (second row of Table |^) and verified that they are compatible with the fitted values. 
The numerical results are consistent with the picture illustrated above, and the value of a 
is significantly different from 1, demonstrating that such an effect cannot be neglected. 



'For this reason, this effect has no consequences for the results presented in this paper. 
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(1.86 X 10^ 


(2.82 X 10^ 


0.66 + 5 


Kh = 0.125 


(2.86 X 10^ 
(2.81 X 10^ 


(2.79 X 10^ 
(2.85 X 10^ 


(1.99 +;^) X 10^ 


(2.83 +5^) X 10^ 


0.71 + I 


Kh = 0.129 


(2.83 +3^) X 10^ 
(2.76 X 10^ 


(2.73 +29) X 
(2.80 + ) X 10^ 


(2.17 +'°) X 10^ 


(2.78 X 10^ 


0.78 + ^ 


Kh = 0.133 


(2.77 +32) X 10^ 
(2.68 +23) X 


(2.62 +^^) X 10^ 
(2.71 +^5) X 10^ 


(2.29 +^^) X 10^ 


(2.69 +3q) X 10^ 


0.85 + I 



Momentum p = ( — ,—, 0) 





72 


72 


72 




a 


Kh = 0.121 


(1.94 +|) X 10^ 
(1.90 +3^) X 10^ 


(1.93 X 10" 
(1.97 +3^) X 10" 


(1.09 X 10" 


(1.94 X 10" 


0.56 


K;, = 0.125 


(1.93 +33) X 10^ 
(1.89 X 10^ 


(1.91 +33) X 10" 
(1.96 +3^) X 10" 


(1.16 X 10" 


(1.92 +33) X 10" 


0.61 

— 


Kh = 0.129 


(1.93 X 10^ 
(1.87 +3') X 10^ 


(1.88 +33) X 10" 
(1.95 +3^) X 10" 


(1.25 X 10" 


(1.91 +3°) X 10" 


0.66 


K,, = 0.133 


(1.89 X 10^ 
(1.80 X 10" 


(1.79 +32) X 10" 
(1.88 X 10" 

^ —3d' 


(1.34 X 10" 


(1.84 +|) X 10" 


0.73 



Table 9: Estimates of a and Zl from the fitted values of Zm, Ze, Zp. In the second row cor- 
responding to each Kh, the fitted Zm and Ze are compared with the estimates using Eq. p^j, 
and the measured values of Z"^ and a. 



We conclude this appendix with a discussion of the SL correlators, whose spin structure is 
again different from that of local correlators. This feature must be taken into account in the 
analysis of three-point correlators when the inserted current operator is local. Following the 
reasoning above, we find: 



Momentum 


Kh = 0.121 


Kh = 0.125 


Kh = 0.129 


Kh = 0.133 


p = 

p = (t.0'0) 


(3.8 +3) X 10-3 
(3.9 +J) X 10-3 
(3.8 X 10-3 


(3.6 +3) X 10-3 
(3.7 X 10-3 
(3.7 +J) X 10-3 


(3.5 +3) X 10-3 
(3.6 +J) X 10-3 
(3.5 X 10-3 


(3.2 X 10"' 
(3.3 +J) X 10-3 
(3.3 +5 X 10-3 



Table 10: Values of Zi = (ZiZ^)/ J , for the four heavy masses and kh = K12 = 0.14144. 



I 4£/ 4l!/ 

(1 + 7o) + a(l - 7o) 



p-7 .(34) 



As above, it is possible to measure ZsZi averaging Ze and Zm- By doing so, we have extracted 
Zi for three different values of the momentum, and the results are shown in Table |l^. The 
evidence that Zi is independent of p is a further check of the validity of our interpretation. 
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Figure 1: Effective masses for the A and E haryons. We show typical plots, corresponding 
to Kh = 0.129 and kh — K12 = 0.14144. The straight lines are our best fits, which agree for 
SL and SS correlators starting from imin = 17 onwards. 
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Figure 2: Masses and x^/dof obtained from a sliding window analysis for the A and the E 
baryon correlators. 




Figure 3: The chiral behaviour of the A and E masses. The boxes denote data at our three 
light-quark masses; the crosses denote the extrapolation of our results to the chiral limit and 
the strange- quark mass. 




Figure 4: Extrapolation of the A and S baryon masses. In the figures both the linear extrap- 
olation and that obtained taking into account 0{1/Mp) correction are shown. The diamonds 
and crosses correspond to the extrapolations to the charm mass and beauty mass respectively. 
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Figure 5: M\ — Mp splitting: comparison of the values obtained by different groups and using 
different fermion actions. The estimates are also compared with the experimental numbers. 
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Figure 6: M\ — Mp splitting as obtained adopting the two procedures A and B. The solid lines 
correspond to the linear extrapolation in the inverse pseudoscalar mass and the dotted line is 
the same extrapolation modelled with a quadratic dependence. In the plot, the light k, values 
are also indicated. The results, which are consistent between both methods, are conpared with 
the experimental values. The vertical dotted lines indicate 1/Md and l/Mp. 
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Figure 7: Ms — M\ as a function of 1/Mp. The linear and the quadratic extrapolation are 
shown. The vertical dotted lines indicate l/Mc and 1/Mb, rspectively. 
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Figure 8: Ms* — My. splitting computed from the mass difference, together with linear and 
quadratic eoctrapolations. The vertical dotted lines indicate l/M^ and 1/Mb- 
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